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A  dynamic  model  for  a  pouch  type  Li-polymer  battery  based  on  electrochemical  and  thermal  principles  is 
developed  to  analyze  static  and  dynamic  performances  of  a  single  cell.  The  model  for  the  single  cell  is 
a  quasi-three-dimensional,  constructed  by  connecting  one-dimensional  models  for  micro  cells  using 
current  collectors.  The  developed  model  can  represent  distributions  of  temperature,  potentials,  and 
current  flows  along  with  distribution  of  lithium  ions  through  the  plane.  The  model  is  coded  using 
MATLAB  and  validated  against  a  LiMn204/Carbon  pouch  type  power  cell.  The  static  analysis  includes 
responses  of  the  terminal  voltage  at  different  current  rates  derived  from  the  distribution  of  over¬ 
potentials  in  the  micro  cell  as  function  of  state  of  charge  (SOC)  as  well  as  distribution  of  potentials  and 
current  flows  in  the  single  cell.  Ion  concentration  in  the  electrodes  and  electrolyte  of  the  micro  cell  are 
analyzed.  The  dynamic  analysis  includes  voltage  and  temperature  responses  during  charging  and  dis¬ 
charging  processes.  The  results  demonstrate  effects  of  operation  conditions  on  key  variables  of  the  cell 
performance  that  includes  distribution  of  ions  in  the  electrodes  and  electrolyte  in  micro  cells  as  well  as 
distribution  of  heat  generation  in  single  cell  level  during  a  charging  and  discharging  process. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Energy  efficiency  of  future  power  systems  can  be  substantially 
improved  using  energy  storage  that  can  capture  energy  and 
retrieve  energy  that  is  otherwise  dissipated.  Batteries  are  the  most 
preferred  energy  storage  because  of  high  columbic  efficiency,  high 
energy  and  power  density.  Among  different  batteries  technologies, 
Li-polymer  battery  is  mostly  preferred  because  of  its  highest  power 
and  energy  density,  and  its  drastic  price  drop  that  was  triggered  by 
the  rapidly  growing  electronic  market.  The  mass  production  of  cells 
for  the  consumer  electronics  ensures  a  high  production  volume 
with  increased  quality.  On  the  other  hand,  new  applications  such  as 
electric  and  hybrid  vehicles  require  increased  power  and  energy 
density  that  is  accomplished  by  employing  a  large  size  of  active 


*  Corresponding  author.  Tel.:  +1  334  844  3329;  fax:  +1  334  844  3307. 

E-mail  addresses:  mzx0001@auburn.edu,  xiaomeng996@sohu.com  (M.  Xiao), 
choe@auburn.edu  (S.-Y.  Choe). 

0378-7753 /$  —  see  front  matter  ©  2012  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.org/10.1016/jjpowsour.2012.05.103 


areas  by  means  of  wrapping  or  folding  electrodes,  electrolytes,  and 
separators  together.  Dependent  upon  the  way  of  packaging,  the 
cells  can  be  manufactured  into  three  different  types,  cylindrical, 
prismatic,  or  pouch  type.  Compared  with  the  cells  used  for 
consumer  electronics,  these  high  power  cells  can  be  quickly 
degraded  and  performances  might  drop  fast  and  becomes  unstable 
because  of  high  ion  transport  and  gradient  of  ion  concentrations 
and  the  associated  high  heat  generation  in  conjunction  with 
varying  environmental  operating  conditions.  The  working  mecha¬ 
nism  of  a  cell  is  very  complex  and  hard  to  understand  in  details, 
which  presents  one  of  technical  barriers  that  prevents  from  optimal 
design  of  a  cell  and  systems. 

Terminal  behaviorof  a  battery  cell  can  be  described  in  three 
different  ways,  using  empirical  equations,  equivalent  circuit 
components,  or  electrochemical-thermal  principles.  The  first  two 
approaches  are  able  to  describe  terminal  behaviors,  but  ignore 
detailed  physical  effects  that  include  spacious  distribution  of 
charging  and  discharging  ionic  and  electronic  current  and  the 
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Nomenclature 

W 

Work  (J) 

a 

Specific  surface  area  for  electrode  (cm-1) 

ti 

Initial  transference  number 

c 

Concentration  (mol  L-1) 

Greek 

symbols 

D 

Diffusion  coefficient  (cm2  s-1) 

e 

Porosity  of  a  porous  medium 

E 

Energy  storage  (J) 

</> 

Potential  in  a  phase  (V) 

F 

Faraday  constant  (96,487  C  mol-1) 

V 

Surface  overpotential  (V) 

e 

Energy  storage  density  (J  cm-1) 

K 

Ionic  conductivity  of  electrolyte  (S  cm"1) 

*0 

Reference  exchange  current  density  (A  cm'2) 

Kd 

Diffusion  conductivity  (A  cm-1) 

j 

Transfer  current  density(A  cm-2) 

P 

Density  (g  cm-3) 

SOC 

n 

State  of  charge 

Amount  of  active  electrode  material 

<7 

Conductivity  of  solid  active  material  (S  cm-1) 

R 

Universal  gas  constant  (8.3143  J  mol-1  K-1) 

Subscripts 

T 

Cell  temperature  (K) 

a 

Anode  reaction 

U 

Standard  potential  affected  by  solid  concentration  (V) 

c 

Cathode  reaction 

X 

Stoichiometric  number  in  negative  electrode  or 

e 

Electrolyte  phase 

coordinate 

eff 

Effective 

y 

Stoichiometric  number  in  positive  electrode  or 
coordinate 

s 

Solid  phase 

associated  temperature  rise  at  different  state  of  charge  (SOC).  In 
addition,  operating  conditions  like  cycling,  ambient  temperature, 
and  storage  times  cannot  be  considered.  By  contrast,  the  model 
based  on  electrochemical-thermal  principles  can  better  represent 
those  major  effects,  but  parameters  necessary  for  the  model  are 
difficult  to  characterize  and  obtain.  The  electrochemical  model 
known  as  a  thin-film  model  was  proposed  by  M.  Doyle,  T.  Fuller, 
and  J.  Newman  [1].  The  model  was  extended  by  reflecting  other 
cathode  materials  like  UC0O2,  LiMn204  [2]  and  LiFeP04  [3].  Later, 
performance  of  cells  for  a  mixture  of  two  different  active  materials 
in  electrodes,  LiyNio.80Coo.15Alo.05O2  and  LiyMn204  was  studied  [4]. 
In  addition,  the  model  was  used  to  investigate  effects  of  material 
properties  on  performance,  which  include  transference  number  [5], 
activity  coefficient  in  electrolyte  [6],  diffusion  coefficient  in  LiPF6 
electrolyte  [7],  and  particle  size  distributions  on  the  discharge 
capacity  [8].  Furthermore,  the  model  was  improved  by  adding  the 
energy  equation  [9]  and  the  efficiency  of  batteries  during  cycling 
was  investigated  [10].  Effects  of  side  reactions  were  also  studied 
[11].  Other  studies  include  relaxation  phenomena  of  a  lithium  ion 
battery  [12],  current  limitation  on  pulse  operation  [13],  degradation 
of  the  battery  in  a  pocket  computer  [14],  calculation  of  heat  transfer 
[15],  and  optimal  design  of  a  battery  system  used  for  electric 
vehicles  [16].  However,  the  research  work  above  has  not  extended 
to  a  high  power  pouch  cell  type  lithium  Polymer  battery  that  is 
currently  widely  accepted  for  applications  to  electric  vehicle  (EV) 
and  hybrid  electric  vehicle  (HEV).  Therefore,  we  developed  a  quasi- 
three-dimensional  dynamic  model  based  on  the  principles  of 
electrochemical  kinetics,  mass  balances,  charge  conservations  and 
energy  balance.  The  model  is  capable  of  representing  ion  transport, 
gradients  of  ion  concentrations  and  two-dimensional  distribution 
of  potentials,  current  density  and  temperature  as  well  as  terminal 
voltage  and  current  as  a  function  of  time. 

2.  Setup  for  development  of  a  model  for  a  single  cell 

A  micro  cell  and  a  pouch  type  of  a  single  cell  are  shown  in  Fig.  1, 
respectively.  The  single  cell  is  assumed  to  be  made  of  micro  cells  that 
are  connected  in  parallel  by  current  collectors  of  both  electrodes. 

Each  of  the  micro  cells  is  assumed  to  be  a  sandwiched  model  that 
is  composed  of  a  negative  electrode,  a  separator,  and  a  positive 
electrode.  The  materials  are  porous  and  equally  mixed  with  elec¬ 
trolyte  materials.  In  addition,  particles  dispersed  at  the  cathode  and 


anode  electrodes  have  a  same  radius,  rs,  and  are  in  contact  with 
neighboring  particles.  The  working  mechanism  of  the  micro  cell  is 
described  using  five  variables  that  present  concentrations  in  elec¬ 
trodes  and  electrolyte,  potentials  in  solid  and  electrolyte,  and  over¬ 
potentials  under  isothermal  condition.  Since  the  conductivity  of 
current  collectors  is  relatively  high,  it  is  assumed  that  no  lateral 
current  flows  from  one  micro  cell  to  another,  so  all  lateral  effects  can 
be  ignored.  In  addition,  the  current  collectors  on  each  layer  have  the 
same  potential,  so  that  only  two  current  collectors  are  considered  for 
the  model.  The  micro  cell  becomes  a  one-dimensional  model.  The 
tab  current  of  the  two  collectors  in  the  model  are  obtained  by 
dividing  total  tab  current  by  the  number  of  layers. 

The  number  of  micro  cells  is  equal  to  that  of  the  grids,  given  by 
meshing  the  single  cell.  The  number  of  grids  in  the  following 
calculation  is  1900,  resulting  from  38  divisions  in  the  x-direction 
(x  =  38)  and  50  divisions  in  the  y-direction  (3/  =  50).  The  area  of 
each  individual  micro  cell  may  be  obtained  by  dividing  the  cell  area 
by  the  number  of  grids. 

When  ions  and  electrons  are  transported  though  electrolyte  and 
external  circuit,  chemical  reactions  will  take  place  at  the  surface  of 
the  electrode  particles  that  contact  electrolyte.  After  the  reactions 
have  been  completed,  the  free  ions  are  located  in  electrodes  and 
diffuse  until  intercalated.  The  electrodes  are  approximated  with 
spheres  and  it  is  assumed  that  the  lithium  ions  diffuse  only  in  the 
radial  direction,  driven  by  the  gradient  of  the  Li-ion  concentration. 

A  schematic  diagram  for  the  model  of  a  single  cell  is  shown  in 
Fig.  2,  where  initial  conditions,  input  variables  and  parameters  for 
the  model,  and  outputs  are  included.  The  parameters  used  for  the 
model  are  appended.  The  input  variables  are  electric  loads  that  can 
be  bi-directional  current  and  voltage  as  a  function  of  time,  while 
the  outputs  are  responses  of  the  loads  in  terminal  voltage  or 
current,  SOC,  and  temperature  distributions  as  well  as  other 
internal  variables.  In  addition,  the  diagram  shows  how  governing 
equations  of  the  micro  cells  are  connected  each  other  and  with 
those  in  the  single  cell. 

3.  Governing  equations 

3.1.  Overpotentials 

Current  generated  by  chemical  reactions  are  described  by  the 
Butler- Volmer  equations,  shown  in  Eq.  (1). 
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Fig.  1.  A  micro  cell  and  a  pouch  type  single  cell. 
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where/1  is  the  current  density,  i'o  is  the  exchange  current  density,  as 
is  the  specific  interfacial  surface  area,  0S  is  the  solid  potential,  0e  is 
the  electrolyte  potential,  U  is  the  equilibrium  potential,  T  is  the 
absolute  temperature,  n  is  the  number  of  electrons  involved  in  the 
electrode  reaction,  F  is  the  Faraday  constant  and  R  is  the  universal  gas 
constant,  a  is  the  symmetry  factor  that  is  dimensionless  and  pres¬ 
ents  a  ratio  between  oxidation  and  reduction,  tjsei  is  the  additional 
overpotential  caused  by  the  Solid  Electrolyte  Interphase  (SEI)  layer. 

The  exchange  current  density,  io,  is  the  intrinsic  rates  of  elec¬ 
trons  and  ions  transfer  at  the  interface  between  electrolyte  and 
electrodes.  It  can  be  expressed  with  concentrations  and  symmetric 
factors,  cta  and  ac,  that  describe  bi-directions  of  one  reversible  semi¬ 
reaction. 


jLi  =  °s 


•'of 


exp 


Tta  -n-F 
RT 


(v-vsm)  -exp 


ac-n- 

“RT 


*o 


max 


(2) 


where  cs  and  ce  refer  to  the  lithium  ion  concentration  in  solid  and 
electrolyte  and  cSimax,  is  the  theoretical  maximum  ion  concentra¬ 
tion  in  solid.  aa  and  ac  is  the  oxidation  and  reduction  symmetry 
factor. 


3.2.  Charge  conservations ,  ions  transport,  and  mass  balance 

The  charges  produced  in  oxidation  processes  should  be  the 
same  as  those  consumed  in  the  reduction  process.  This  relationship 
is  described  using  Ohm’s  law,  expressed  as  a  function  of  current  and 
potential  gradients.  Ion  transport  in  the  electrolyte  is  governed  by 
the  Nernst-Planck  equation,  where  migration  and  diffusion  terms 
are  considered,  as  shown  in  the  following  equation. 

^  (/■!/)  +!/^Z(lnce)) +jU  =  0  (3) 


Initial  conditions: 
•Initial  SOC 
•Initial  temperature 
distribution 


Input  variables: 
•Load  profile 
•Ambient 
temperature 


Parameters: 

•Battery  geometry 
•Maximum  capacity 
•Concentration 
•Activity  coefficient 
•Diffusion 
coefficient 
•Change  of  enthalpy 
•Conductivity 

etc. 
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Fig.  2.  Scheme  of  the  model. 
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where/1  is  the  current  density  in  the  electrolyte  or  solid  phase,  0e  is 
the  electrolyte  potential,  ce  is  the  concentration  of  the  electrolyte, 
Keff  is  the  effective  electrolyte  conductivity,  while  /coff  is  the 
concentration  driven  diffusion  conductivity. 

As  for  charge  transport  in  solid,  no  diffusion  term  is  necessary 
for  electron  transport  in  solid  because  it  is  solely  governed  by 
Ohm’s  law. 


0 

0X 


where  creff  is  the  solid  phase  electronic  conductivity. 


(4) 


3.3.  Material  balance 


Charges  produced  or  consumed  by  chemical  reactions  taking  place 
in  the  electrodes  are  conserved.  The  chemical  process  in  the  entire  cell 
can  be  described  using  a  material  balance  principle.  Hence,  the  balance 
of  materials  in  electrolyte  is  affected  by  the  gradient  of  ion  concentra¬ 
tions.  Since  the  electrode  pores  are  filled  with  electrolyte,  the  balance 
should  consider  the  porosity  of  the  material.  Therefore,  the  diffusion 
coefficient  is  redefined  by  considering  the  porosity  as  follows: 

of  =  De-e<?  (5) 

where  the  diffusion  coefficient,  De,  is  constant.  The  material  balance 
equation  for  the  electrolyte  considering  the  porosity  is  as  follows: 


d(eeCe) 

dt 


+ 


1 


0 


(6) 


where  ce  is  the  concentration  of  lithium  ion  in  electrolyte,  F  is  the 
Faraday  constant,  t+  is  the  initial  transference  number,  ee  is  the 
porosity,  and  jLl  is  the  current  density  [17]. 

For  electrodes, 


0CS  _  Ds  0  /  2^s\ 
dt  r2  dr  \  dr ) 


(7) 


where  cs  is  the  lithium  ion  concentration  in  electrode  and  r  is  the 
radius  of  an  electrode  particle. 

When  reaction  rates  change,  the  current  density  varies  accord¬ 
ingly.  Consequently,  the  concentration  is  changed. 

The  reaction  rate,  jLl,  determined  by  the  Butler-Volmer  equa¬ 
tion  is  the  most  important  variable  that  affects  not  only  the  charge 
conservation  equations  but  also  the  mass  transport  equation. 


3.4.  Energy  conservation 


Temperature  in  operating  cells  affects  the  performance  and  degra¬ 
dation  of  materials,  which  is  described  using  the  energy  equation: 


dT 

dt 


0 

0X 


where  p,  Cp,  k  and  Q.  are  the  density,  heat  capacity,  thermal 
conductivity,  and  heat-generation  rate  per  unit  volume, 
respectively. 

The  heat  source  term  below  in  Eq.  (9),  Q,  is  generally  used  to 
calculate  the  temperature  [18]. 


Q 


^ocv  - 


T- 


dUocv\ 
dT  ) 


(9) 


currents  that  are  caused  by  the  non-uniformity  of  electrochemical 
energy  state  gradients,  which  is  called  the  heat  of  mixing.  There¬ 
fore,  the  heat  source  term,  Q,  is  divided  into  three  terms,  a  revers¬ 
ible  heat  in  micro  cell  (qrev),  an  irreversible  heat  in  micro  cell  (girr), 
and  the  Joule  heating  caused  by  the  current  collectors  (qcc)  in 
a  single  cell  per  unit  volume. 


Q  —  J  QrevdV  +  J  Qirr.dV  +  J  Qcc.dV 


V 


V 


V 


(10) 


The  reversible  heat  is  a  result  of  the  change  of  entropy  during 
a  chemical  reaction  (FAS).  The  rate  of  the  heat  generated  per 
volume  can  be  expressed  as  follows: 

r\F 

qi-ev  =jUT—  (11) 

For  calculation  of  irreversible  heat  generation,  a  following 
modification  is  suggested.  When  a  battery  is  charged,  electric 
energy  supplied  to  the  cell  at  the  terminal  is  equal  to  a  sum  of 
electrochemical  energy  stored  in  the  battery  and  irreversible  heat 
dissipated.  When  discharged,  the  electric  energy  obtained  at  the 
terminal  is  the  difference  between  the  electrochemical  energy  and 
the  irreversible  heat.  When  the  charging  current  is  very  low,  it  can 
be  assumed  that  the  energy  dissipated  is  negligible.  Then,  the 
electrochemical  energy  can  be  expressed  as  follows: 


soc 

—  Qmax  J  UoQydSOC 

SOCo 

SOC 

=  Qmax  J  (u+qu  -uequ)dsoc 

SOCo 


(12) 


Since  there  are  always  gradients  of  lithium  ion  concentrations, 
equilibrium  potentials  and  current  density,  total  amount  of  the 
electrochemical  energy  can  be  obtained  by  integrating  the  energy 
in  a  small  control  volume.  On  the  other  hand,  change  of  SOC  is 
proportional  to  change  of  concentration  in  a  volume,  so  the  SOC  in 
the  equation  above  can  be  substituted  by  Cs  as  follows: 


Chemical  =  -^s+F- J  Cs+  ■  UeqU  ■  dV  -  eS-  ■  F  ■  J  CS--UequdV  (13) 


V 


V 


where  es  is  the  active  material  volume  fraction,  F  is  the  Faraday 
constant,  cs  is  the  solid  phase  lithium  ion  concentration  and  Uequ  is 
the  equilibrium  potential. 

The  irreversible  heat  generation  rate  for  a  micro  cell  is  the 
difference  between  the  electrical  power  supplied  to  the  cell 
(— Wnicro  •  Wro)  and  the  power  (Pchem.)  contributing  to  the  chem¬ 
ical  energy  increase  that  is  obtained  by  differentiating  the  elec¬ 
trochemical  energy  with  respect  to  the  time.  When  charging, 
0  <  Pchem.  <  —  ^micro  *  fmicro  and  when  discharging,  Pcheml  <  —  ^micro 
•  Wro<0,  so  that  qirr  always  becomes  positive. 


<?irr  =  ~Y 


micro  1  micro 


-P 


chemical 


(14) 


The  heat  generation  rate  caused  by  the  Joule  heating  in  the 
current  collectors  is  as  follows: 


qcc  =  <t-(V<A cc)2  (15) 

In  addition,  the  heat  transfer  by  convection  is  considered. 

Qconvec.  —  *  (Ps  —  Too )  (16) 


The  equation  indicates  that  no  heat  will  be  generated  when  the 
current  goes  zero.  In  reality,  heat  is  continuously  generated  by  the 


where  h,  Is,  and  Too  denote  the  convective  heat  transfer  coefficient, 
the  surface  temperature,  and  ambient  temperature,  respectively. 
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3.5.  Model  simplifications 


Terminal  voltage  of  a  micro  cell  is  given  as  a  difference  between 
electrode  potentials  and  equilibrium  potential.  The  electrode 
potentials  are  a  function  of  concentration  given  by  the  Nernst 
equation.  The  equilibrium  potential  of  a  cell  is  called  open  circuit 
voltage  (OCV)  when  no  current  flows.  The  OCV  is  the  difference 
between  the  two  standard  potentials  of  the  positive  and  negative 
electrode  and  is  measured  experimentally.  Separation  of  individual 
equilibrium  potentials  for  positive  and  negative  electrodes  from  the 
measured  OCV  are  carried  out  using  an  empirical  equation  for  the 
equilibrium  potential  in  negative  electrode  given  in  [19]. 


U-(x)  =  8.00229+5.0647x— 12.578X1/2  — 8.6322e— 4x_1 

+2.1 765e-5x3/2  -  0.4601 6exp(l  5(0.06  -x)) 
-0.55364exp(— 2.4326(x— 0.92))  (17) 

Currents  produced  during  chemical  reactions  are  a  function  of 
overpotentials  governed  by  the  Butler-Volmer  equation  that 
consists  of  two  exponential  functions.  In  fact,  the  overpotentials 
usually  vary  within  a  linear  range  under  normal  battery  operating 
conditions,  Therefore,  the  exponential  function  can  be  linearized  as 
shown  in  Eq.  (18).  Comparison  between  the  exponential  and  the 
linearized  function  shows  that  the  difference  is  about  0.02-0.04  V 
for  the  overpotential  value  up  to  0.4  V.  Particularly,  the  use  of  this 
linear  equation  in  the  charge  conservation  equation  eases  conver¬ 
gences  of  all  calculations  and  reduces  computational  time. 


.  ($a  +  &c)F 

asl°S — rT^v 


(18) 


As  matter  of  fact,  the  time  constant  of  ion  concentration  is  larger 
than  that  that  for  electrons,  it  is  assumed  that  the  second  term  in 
the  conservation  of  charge  Eq.  (3)  is  regarded  as  constant.  There¬ 
fore,  concentration  and  temperature  at  a  given  time  can  be  used  to 
calculate  /ceff  and  /cof.  Hence,  the  only  unknown  variables  are  the 
phase  potentials,  so  that  the  equations  are  linearized  and  the 
number  of  equations  to  be  solved  simultaneously  is  reduced. 


3.6.  Numerical  calculation 


3.6.1.  Initial  and  boundary  conditions 

Numerical  calculation  requires  initial  and  boundary  conditions 
to  initiate  and  finalize  simulations.  The  three  variables  calculated 
from  dynamic  equations  are  lithium  ion  concentration  in  solid  and 
electrolyte,  and  cell  temperature,  while  others  are  obtained  from 
static  equations  using  these  three  variables. 

The  boundary  condition  for  the  potential  of  an  electrode  is: 

effd^s  _  effd0s 
-  Ox  x=0_  +  dx 

where  the  creff  is  the  conductivity,  0S  is  the  solid  phase  potential,  x  =  0, 
L  represents  the  interface  between  the  electrode  and  current  collector 
at  negative  and  positive  sides,  respectively,  and  x  =  <5_,  5+  represents 
the  interface  between  the  separator  and  electrode  at  negative  and 
positive  sides,  respectively.  For  a  given  current  density,  the  potentials 
outside  of  the  electrolyte  boundary  are  zero  as  follows: 


x=L 
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90s 


Ox 
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x=S 


Ox 


=  0  (19) 
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The  initial  conditions  of  various  variables  are  determined  from 
the  initial  lithium  ion  concentration  in  the  electrodes.  The 
concentration  is  expressed  as  a  function  of  SOC,  as  follows: 
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x=L 


(20) 


cs_  =  stoi0  +  stoi 
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Cs+  —  (stoij  +  (stoi|00  —  Stoij)  X  SOC^  ’Cs+niax 


(21) 


The  concentration  of  lithium  ions  inside  the  particles  is  spherically 
symmetric  with  no  generation  within,  so  the  gradient  of  lithium  ion 
concentration  at  the  inner  boundary  is  zero.  The  outer  boundary  of 
electrode  particles  is  determined  by  the  rate  of  the  reaction  taking  place. 
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r=Rs 


acf 


Since  there  is  no  electrolyte  flow  prior  to  current  flow  through 
the  battery,  the  change  in  electrolyte  concentration  at  the  bound¬ 
aries  is  shown  below  to  equal  zero. 
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(23) 


The  initial  cell  temperature  is  set  to  an  ambient  temperature.  At 
a  low  current,  heat  transfer  between  the  environment  and  the  cell 
is  regarded  as  negligible.  At  the  boundary, 

-fcg  =  ft-(r-Ta)  (24) 

where  k  is  the  thermal  conductivity  of  the  outer  casing,  h  is  the 
overall  convective  heat  transfer  coefficient,  Ta  is  the  ambient 
temperature,  and  T  is  the  temperature  in  the  cell. 


3.6.2.  Numerical  method 

The  micro  cell  is  discretized  into  one  hundred  grids  through  the 
plane.  The  negative  and  positive  electrodes  are  meshed  with  45  and 
34  grids,  respectively,  because  of  different  thickness  of  the  elec¬ 
trodes.  As  a  result,  179  discrete  equations  are  obtained  for  the  Eqs.  (3) 
and  (4)  in  order  to  calculate  potentials.  As  potentials  and  concen¬ 
trations  near  both  interfaces  between  electrode  areas  and  separator 
change  drastically,  the  size  of  grids  at  the  interfaces  should  be  small 
enough  to  minimize  the  numerical  errors.  In  addition,  the  small  grid 
is  generated  at  the  outer  boundary  of  spheres. 

All  of  equations  are  converted  into  discrete  forms,  coded  using 
MATLAB  and  solved  for  the  solid  phase  potentials  at  the  electrodes 
and  the  electrolyte  phase  potential  across  the  cell.  For  calculations, 
a  finite  differential  method  with  an  implicit  scheme  is  employed  to 
get  converged  fast. 

The  resulting  phase  and  equilibrium  potentials  are  used  to  calcu¬ 
late  overpotential  at  each  section  across  the  cell,  which  allows  for 
calculation  the  reaction  rate,  jLl,  using  the  Butler-Volmer  equation. 
The  reaction  rate  will  be  inputted  to  mass  balance  Eqs.  (6)  and  (7)  to 
determine  the  change  in  concentration  at  each  time  step.  Similar  to 
the  method  used  for  the  calculation  of  the  reaction  rate,  an  equation 
for  each  segment  or  grid  of  the  battery  cell  is  derived.  While  the  initial 
values  are  known,  the  lithium  ion  concentration  as  a  function  of  the 
change  in  time  provides  the  variables  for  the  equation. 

In  fact,  electric  field  in  both  of  electrodes  and  electrolyte 
responses  faster  than  that  of  ion  transports  when  charged  and 
discharged.  Therefore,  it  is  assumed  that  the  potentials,  4>s  and  4>e 
are  constant  during  each  time  step  and  calculated  just  in  one 
subroutine.  Conversely,  ion  transport  is  dynamically  calculated  in 
another  subroutine  using  the  mass  balance  equations  for  ion 
concentrations  cs  and  ce.  The  discretized  equations  can  be  numer¬ 
ically  solved  by  implicit  or  explicit  method.  The  major  drawback  of 
the  implicit  method  is  the  difficulty  in  coding  in  compared  with  the 
explicit  one  that  is  less  accurate  and  stable.  The  code  developed  is 
based  on  the  implicit  method.  The  dynamic  equation  above  has 
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a  general  form  dc/dt  =/(c,0),  where  c 


and  (f> 


The 


potential  can  be  solved  from  the  static  Eqs.  (3)  and  (4).  Even  though 
both  explicit  and  implicit  method  assume  that  the  function,  f 
during  two  time  steps  is  constant,  the  explicit  method  explicitly 
calculates  the  functionusing  c  and  f  from  the  previous  time  step, 
while  the  implicit  method  considers  c  and  f  to  be  the  values  of 
potentials  and  concentrations  in  the  following  time  step  and 
calculate  them  by  solving  a  matrix  about  c.  For  calculation  of  the 
potential,  at  the  following  time  step,  change  of  concentration  for 
the  time  step,  Ac,  instead  of  c  is  calculated  at  a  Co  at  the  previous 
time  point,  which  results  in  the  potential  f  as  follows; 


potential  is  relatively  small  and  varies  from  x  =  0  to  x  =  L.  In  fact,  the 
activation  overpotential  is  the  difference  between  the  potential  in 
solid,  electrolyte,  and  equilibrium  potentials,  as  shown  in  Eq.  (17). 

Fig.  4  (c)  shows  the  cathode  and  corresponding  equilibrium 
potential  and  overpotential.  The  overpotential  at  the  cathode 
always  has  an  opposite  sign  when  compared  to  the  anode,  and  in 
each  electrode  it  can  be  either  positive  or  negative  depending  on 
whether  the  cell  is  charging  or  discharging.  In  this  case,  the  cell  is 
obviously  discharging  because  the  overpotential  at  the  anode  is 
positive,  which  means  the  current  is  flowing  from  the  negative 
electrode  to  the  electrolyte  then  to  the  positive  electrode.  Fig.  4  (d) 
shows  the  corresponding  overpotential.  Fig.  5 


4>  =  4>o  +  A-Ace  +  BAcse  (25) 

where  the  f o  is  the  potential  at  previous  time  step,  and  A  and  B  are 
coefficient  matrices  that  include  the  linearized  relationship 
between  cj)  and  Ac. 

4.  Static  and  dynamic  analysis  and  validation  of  a  micro  cell 

4.1.  Potential  and  overpotential 

Standard  potentials  are  determined  by  the  type  of  chemical 
compounds  used  for  the  electrodes  and  are  a  function  of  the  stoi¬ 
chiometric  numbers  of  the  chemical  compound  and  SOC.  The 
terminal  voltage  for  a  given  electrode  can  be  calculated  based  on  the 
standard  potentials.  When  a  battery  is  charged  and  discharged,  two 
reversible  half  reactions  occur  that  cause  the  potential  losses 
(Overpotentials).  The  overpotentials  of  a  cell  at  two  different  SOCs 
(50%  and  80%)  are  shown  in  Fig.  3  under  assumption  that  the  OCV  at 
a  specified  SOC  is  constant.  The  concentration  overpotential,  A(J)e,  is 
the  potential  difference  in  the  electrolyte  phase  caused  by  the 
concentration  of  the  reactants.  In  addition,  the  ohmic  and  the 
concentration  overpotential  are  depedent  upon  the  geometry  and 
tend  to  increase  when  the  film  layers  become  thicker.  Conversely, 
the  activation  overpotential  was  not  “cumulative”  through  the  plane 
of  the  battery  film,  so  that  the  value  is  relatively  small.  The  idea  cell 
voltage  in  each  figure  is  constant,  but  varies  as  a  function  of  SOC. 

The  Fig.  4  (a)  shows  the  potential  of  a  micro  cell,  where  the  x 
coordinate  is  a  dimensionless  number  of  meshed  grids  along 
through-the-plane.  The  range  of  the  grid  for  0  <  x  <  45, 46  <  x  <  67, 
68  <  x  <  100  corresponds  to  the  anode,  the  separator,  and  the 
cathode,  respectively.  We  assumed  that  the  potential  of  the  elec¬ 
trolyte  at  x  =  0  is  the  ground  potential,  so  that  the  potential  of  the 
anode  and  the  cathode  is  about  0.1  V  and  3.6  V,  respectively.  The 
zoomed-in  potentials  in  the  anode  are  shown  in  Fig.  4  (b).  The 
difference  between  the  equilibrium  potential,  U,  and  the  solid 


4.2.  Current  distribution 

When  a  battery  is  charged  or  discharged,  chemical  reactions  take 
place  at  the  interface  between  the  electrodes  and  electrolyte.  As 
a  result,  electrons  and  ions  are  separated  at  the  cathode.  The  electrons 
and  ions  are  transported  through  an  external  circuit  and  the  electro¬ 
lyte, respectively,  to  anode  and  combine  again  there  to  complete  the 
total  reaction.  The  current  in  the  electrodes  is  composed  of  electron 
current  and  ion  current  because  of  the  porosity,  while  only  ion  current 
can  go  through  the  separator.  The  gradient  of  ion  currents  is  the  current 
density,  jLl,  whose  slope  indicates  the  direction  of  the  current  (i.e. 
positive  slope  implies  that  the  current  flows  from  solid  to  electrolyte). 

4.3.  Concentration  of  lithium  ions 

Concentration  of  lithium  ions  in  solids  and  electrolyte  is  the  key 
factor  that  affects  physical  variables  and  states  of  the  battery  like  SOC, 
conductivities,  and  exchange  current  density.  Since  the  concentration 
in  the  solid  is  not  uniformly  distributed,  only  the  lithium  ion  concen¬ 
tration  in  the  electrodes  at  the  outer  boundary  is  shown  in  Fig.  6,  which 
determines  the  maximum  discharge  or  charge  rate  at  a  moment. 

Responses  of  a  high  discharging  current  rate  (10C  rate)  on 
concentration  and  current  density  are  shown  in  Fig.  7.  The  current 
density  shows  a  highest  value  in  the  electrode  at  the  point  of 
contact  with  the  separator.  It  decays  with  time  because  the 
migration  of  ions  in  both  the  electrolyte  and  the  electrodes  miti¬ 
gates  the  difference  of  reaction  rates.  It  takes  several  minutes  for 
ion  concentration  to  be  uniformly  distributed. 

In  Fig.  8,  r  =  0  represents  the  center  of  the  electrode  particles, 
while  r  =  1  is  the  outer  boundary  of  the  electrode  particles.  The 
geometry  of  electrodes  along  with  material  properties  like  thickness 
and  porosity  result  in  different  concentrations,  for  example  at  L  =  0 
and  1  =  1.  The  magnitudes  of  concentrations  at  the  boundaries  of  the 
separators  on  the  negative  electrode  are  different  from  that  on  the 
positive  electrode  because  of  the  high  conductivity. 


Fig.  3.  Overpotentials  of  micro  cell  as  a  function  of  discharging  current  at  different  SOC. 
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current  in  micro  cell 


a  discharging  current  is  applied,  the  terminal  voltage  drops  instan¬ 
taneously  because  of  ohmic  resistances  present  in  the  cell.  Then,  the 
voltage  drops  linearly  because  of  the  change  of  the  OCV  induced  by 
the  change  of  SOC.  At  the  end  of  each  discharge,  the  experimental 
shows  that  the  change  of  the  OCV  becomes  large  and  as  a  result  the 
voltage  drops  drastically.  In  addition,  it  should  be  noted  that  all  of  the 
charges  discharged  is  a  function  of  the  current  rate.  For  example, 
when  discharged  at  1C,  93%  of  the  maximum  capacity  is  released, 
while  85%  of  capacity  is  released  when  discharged  at  IOC. 

Generally,  charging  a  cell  is  carried  out  by  two  modes,  a  constant 
current  mode  (CC  mode)  and  a  constant  voltage  mode  (CV  mode). 
Transition  from  one  mode  to  another  is  simply  implemented  by 
changing  the  boundary  condition  of  the  solid  conductive  Eq.  (4). 
The  boundary  conditions  for  the  two  modes  are  as  follows: 


|  £ju  =  o 

I  0s+  —  0s-  =  V micro 


(26) 


Fig.  5.  Current  distribution  of  micro  cell. 

4.4.  Discharging  and  charging  behavior  along  with  validation 

Simulation  and  experiment  results  of  discharging  characteristics 
of  a  pouch  cell  are  compared  in  Fig.  9,  where  the  terminal  voltages  at 
different  current  rates  are  plotted  as  a  function  of  time.  The  pouch 
cell  used  for  the  experiments  has  a  capacity  of  15.7  A  h  and  its 
dimensions  of  the  active  area  are  149.2  mm  by  197  mm.  When 


£jLi  =  o 

£/-'  =  Wro/Af. 


(27) 


The  cell  is  firstly  charged  with  the  CC  mode  until  the  terminal 
voltage  reaches  4.2  V  and  then  a  constant  terminal  voltage  is 
applied  until  the  cell  is  completely  charged.  While  CV  charging,  the 
current  decays  slowly.  Simulation  and  experiment  results  are 
compared  in  Fig.  10,  where  four  different  current  rates  of  20  A,  40  A, 
60  A  and  80  A  are  applied. 
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Fig.  6.  Lithium  ion  concentrations  and  exchange  current  density  in  micro  cell. 
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Concentration  in  electrolyte  Current  density 


Fig.  7.  Unsteady  behaviors  of  concentration  and  current  at  IOC  discharging. 


Comparison  of  discharging  and  charging  characteristics  at  the 
different  C  rates  shows  that  the  model  developed  can  predicts  the 
performances. 


5.  Modeling  of  a  pouch  type  single  cell 


Micro  cells  are  connected  in  parallel  by  the  current  collectors 
that  provide  a  current  pathway  to  the  terminal  tabs.  For  each  of  the 
current  collectors,  the  potential  field  is  described  using  Ohm’s  law: 
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(28) 


where  er_  and  a+  are  the  conductivity  of  aluminum  and  copper, 
respectively,  0CC-  and  0CC+  are  the  potentials  on  two  current 


collectors,  Jm icro  is  the  current  for  the  micro  cells  with  units  of 
A  cnrT2,  and  /cc_  and  lcc+  are  the  thickness  of  the  current  collectors. 

Since  the  cell  is  made  of  identical  multiple  layers  of  current 
collectors,  it  is  assumed  that  the  collectors  have  the  same  boundary 
conditions  as  the  load  profile  at  the  terminal  tabs.  Wro  that  is  given 
by  derivatives  of  the  potential  field  0CC  on  both  current  collectors,  as 
described  in  Eq.  (28).  On  the  other  hand,  the  relationship  between 
Wro  and  the  voltage  of  the  corresponding  micro  cell, 
VWro  =  W+  -  0cc-.  should  satisfy  the  nonlinear  dynamic  micro  cell 
equations,  so  that  all  of  the  micro  cells  can  be  coupled  to  each  other. 

The  simulation  results  reveal  that  the  /micro  is  asymmetrically 
distributed  because  of  the  different  conductivity  of  the  two  current 
collectors,  even  though  the  terminal  tabs  are  symmetrically  located. 
Responses  are  calculated  at  1  s  after  a  100  A  discharging  current  is 
applied  where  the  initial  temperature  is  300  K  and  the  initial  SOC  is  50%. 

The  potential  distributions  on  both  current  collectors  are  shown 
in  Fig.  11  and  Fig.  12.  At  the  positive  current  collector,  the  potential 
near  the  current  tab  decreases  when  the  battery  is  discharging  and 
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Fig.  8.  Response  of  ion  concentration  in  electrodes. 
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Terminal  voltage  comparison  @  1C/2C/5C  discharge 


Fig.  9.  Discharge  characteristics  at  different  current  rates  (Ambient 
Temperature  =  294  K  and  SOC  =  100%). 


Terminal  voltage  comparison  @  20A/40A/60A/80A  charge 


Current  comparison  @  20A/40A/60A/80A  charge 


Fig.  10.  Charge  characteristics  at  different  rates  in  CC  and  CV  mode. 

increases  when  charging.  For  each  grid,  the  potential  difference 
between  the  two  current  collectors  is  the  potential  of  micro  cell,  so 
the  micro  cells  near  the  two  current  tabs  show  highest  overpotential 
than  other  locations  because  of  a  high  reaction  rate  near  the  tabs. 

Due  to  the  non-uniform  potential  distribution  on  both  current 
collectors,  the  voltages  of  micro  cells  are  not  identical  and  conse¬ 
quently  the  reaction  rates  of  micro  cells  are  different  from  each 
other.  The  reaction  rate  of  micro  cells  determined  by  the  current 
flow  through  the  micro  cells  per  unit  area  with  an  X-Y  cross- 
section  is  shown  in  Fig.  13.  The  reaction  rate  near  terminal  tab  is 


Potential  on  positive  current  collector 


Fig.  11.  Potential  distributions  on  positive  current  collector  when  discharged. 


Potential  on  negative  current  collector 


Fig.  12.  Potential  distributions  on  negative  current  collector  when  discharged. 
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Fig.  13.  Initial  reaction  rates  of  different  micro  cells. 


Terminal  voltage  comparison  for  5C  cycle 
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Fig.  15.  Temperature  distribution  of  the  cell. 


higher  than  that  in  other  area  at  the  beginning.  As  a  result,  the  SOC 
of  the  region  is  low  for  a  period  of  time.  Since  the  low  SOC  leads  to 
a  low  equilibrium  potential,  the  overpotential  decreases  if  the 
terminal  voltage  remains  the  same  as  before  and  consequently  the 
reaction  rate  get  decreased.  The  results  show  that  the  reaction  rate 
near  the  positive  tab  is  the  highest,  20%  greater  than  the  lowest  grid 
at  the  given  operating  condition. 

5.2.  Dynamic  analysis 

In  order  to  assess  dynamic  responses  of  the  model,  an  operating 
condition  for  a  single  cycle  is  defined  with  initial  SOC  of  41%  and 
discharge  current  rate  is  5C  for  200  s,  rest  200  s,  charge  at  5C 
current  for  200  s  then  rest  200  s.  Comparison  of  simulation  and 
experiment  voltage  response  is  shown  in  Fig.  14. 

In  addition,  temperature  distribution  between  experiment  and 
simulation  is  compared  in  Fig.  15,  where  the  temperature  of  the  cell 
is  captured  using  IR  camera.  The  temperature  distribution  in  the 
cell  is  asymmetric  because  of  the  different  conductivities  of  two 
current  collectors.  Temperature  rise  around  the  positive  tab  is 
higher  than  that  that  in  the  negative  tab.  In  addition,  more  heat  is 
generated  during  discharge  process  because  of  exothermic  reac¬ 
tions.  The  temperature  increase  during  discharging  is  always 
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greater  than  that  of  charging  at  a  given  current  rate.  The  calculation 
for  individual  heat  source  terms  shows  a  distribution  of  Joule 
heating  in  electrolyte  (38%),  heat  of  mixing  (24%),  heat  caused  by 
activation  losses  (16%),  change  of  entropy  (10%),  Joule  heating  at 
SEI(8%),  Joule  heating  at  the  current  collectors  (4%)  and  others. 

In  addition,  temperature  of  three  locations  of  the  cell  are 
measured  and  compared  with  that  of  the  simulation  as  shown  in 
Fig.  16.  The  simulation  results  follow  the  experiment  results  with 
some  time  delay  and  errors  that  are  possibly  affected  by  the  heat 
transfer  between  the  cell  and  air  flow  in  a  thermal  chamber. 

6.  Conclusion 

We  developed  a  high  resolution  quasi-three  dimensional  model 
for  a  pouch  type  of  lithium  polymer  cell  with  a  chemistry  of 
LiMn204/Carbon  for  electrodes.  The  model  is  constructed  using 
a  one  dimensional  micro  cell  based  on  electrochemical  and  thermal 
principles  and  experimentally  validated  against  static  and  dynamic 
performance  of  a  cell. 

Analysis  shows  following  interesting  findings: 

1.  Due  to  different  conductivities  of  the  two  current  collectors, 
asymmetrical  distribution  of  the  current  and  voltage  field  is 
observed.  The  asymmetric  current  field  causes  non-uniform 
heat  generation  rates  on  the  current  collectors.  Particularly, 
more  heat  is  generated  on  the  current  collector  located  near  the 
terminal  tabs.  Similarly,  the  asymmetric  voltage  field  affects 
the  reaction  rates  of  micro  cells,  so  heat  generation  rates  and 
the  associated  DOD  are  different  dependent  upon  location. 

2.  Contribution  of  individual  heat  source  terms  are  analyzed, 
where  the  heat  caused  by  Joule  heating  in  electrolyte  and  heat 
of  mixing  are  the  major  parts  at  the  5  C  rate. 

The  model  developed  is  capable  of  representing  not  only  current 
and  voltage  at  the  terminal  tabs  of  a  single  cell  but  also  distribution 
of  ion  concentrations,  potentials,  internal  currents,  and  tempera¬ 
ture  when  a  battery  is  charged  or  discharged.  Future  work  includes 
measurement  of  the  individual  hear  source  terms  and  validation  of 
the  model  for  other  ambient  temperature  ranges. 
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Appendix,  parameters 


Parameter 

Negative  electrode 

Separator 

Positive  electrode 

unit 

Design  specifications  (geometry  and 

Thickness,  5 

50*10 1 2 3  4 5 

25.4*10~4 

36.4*10~4 

cm 

volume  fractions) 

Particle  radius,  Rs 

1 

O 

* 

1 

o 

* 

cm 

Active  material  volume  fraction,  es 

0.58 

0.5 

Polymer  phase  volume  fraction,  sp 

0.048 

0.5 

0.11 

Conductive  filler  volume  fraction,  ef 

0.04 

0.06 

Porosity,  ee 

0.332 

0.5 

0.33 

Solid  and  electrolyte  phase  Li+ 

Maximum  solid  phase  concentration, 

42.2*1  O'3 

74.9*10~3 

mol  cur3 

concentration 

Cs.max 

Stoichiometry  at  0%  SOC:  x0%,  yo% 

0.123 

0.939 

Stoichiometry  at  100%  SOC:  Xi0o%,  y ioo% 

0.685 

0.434 

Average  electrolyte  concentration,  ce 

1.2*10~3 

1.2*10~3 

1.2*1  O'3 

mol  citT3 

Kinetic  and  transport  properties 

Exchange  current  density  coefficient,  k[0 

13.2 

6.79 

A  cm  2 

Charge-transfer  coefficient,  <xa,  ac 

0.5,  0.5 

0.5,  0.5 

Solid  phase  Li  diffusion  coefficient,  Ds 

1.0*10'n 

1.85*10'n 

2  -1 

cm  s 

Solid  phase  conductivity,  a 

1 

0.1 

S  cm  1 

Electrolyte  phase  Li+  diffusion  coefficient,  De 

2.6*10~6 7 8 9 

2.6*1 0~6 

2.6*1 0~6 

2  -1 
cm  s 

Bruggeman’s  porosity  exponent,  p 

1.5 

1.5 

1.5 

Electrolyte  phase  ionic  conductivity,  k 

7.9ce 

7.9ce 

S  cm  1 

exp(-13472Ce'4) 

exp(-13472Ce-4) 

Li+  transference  number,  t° 

0.363 

0.363 

0.363 
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